3.550 \(\int \frac{\sqrt{\cos (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{3/2}} \, dx\)
Optimal. Leaf size=176 \[ -\frac{(7 A-3 B) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} \tanh ^{-1}\left (\frac{\sqrt{a} \sin (c+d x) \sqrt{\sec (c+d x)}}{\sqrt{2} \sqrt{a \sec (c+d x)+a}}\right )}{2 \sqrt{2} a^{3/2} d}+\frac{(5 A-B) \sin (c+d x)}{2 a d \sqrt{\cos (c+d x)} \sqrt{a \sec (c+d x)+a}}-\frac{(A-B) \sin (c+d x)}{2 d \sqrt{\cos (c+d x)} (a \sec (c+d x)+a)^{3/2}} \]
[Out]
-((7*A - 3*B)*ArcTanh[(Sqrt[a]*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(Sqrt[2]*Sqrt[a + a*Sec[c + d*x]])]*Sqrt[Cos[c
+ d*x]]*Sqrt[Sec[c + d*x]])/(2*Sqrt[2]*a^(3/2)*d) - ((A - B)*Sin[c + d*x])/(2*d*Sqrt[Cos[c + d*x]]*(a + a*Sec
[c + d*x])^(3/2)) + ((5*A - B)*Sin[c + d*x])/(2*a*d*Sqrt[Cos[c + d*x]]*Sqrt[a + a*Sec[c + d*x]])
________________________________________________________________________________________
Rubi [A] time = 0.489672, antiderivative size = 176, normalized size of antiderivative = 1.,
number of steps used = 5, number of rules used = 5, integrand size = 35, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.143, Rules used =
{2955, 4020, 4013, 3808, 206} \[ -\frac{(7 A-3 B) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)} \tanh ^{-1}\left (\frac{\sqrt{a} \sin (c+d x) \sqrt{\sec (c+d x)}}{\sqrt{2} \sqrt{a \sec (c+d x)+a}}\right )}{2 \sqrt{2} a^{3/2} d}+\frac{(5 A-B) \sin (c+d x)}{2 a d \sqrt{\cos (c+d x)} \sqrt{a \sec (c+d x)+a}}-\frac{(A-B) \sin (c+d x)}{2 d \sqrt{\cos (c+d x)} (a \sec (c+d x)+a)^{3/2}} \]
Antiderivative was successfully verified.
[In]
Int[(Sqrt[Cos[c + d*x]]*(A + B*Sec[c + d*x]))/(a + a*Sec[c + d*x])^(3/2),x]
[Out]
-((7*A - 3*B)*ArcTanh[(Sqrt[a]*Sqrt[Sec[c + d*x]]*Sin[c + d*x])/(Sqrt[2]*Sqrt[a + a*Sec[c + d*x]])]*Sqrt[Cos[c
+ d*x]]*Sqrt[Sec[c + d*x]])/(2*Sqrt[2]*a^(3/2)*d) - ((A - B)*Sin[c + d*x])/(2*d*Sqrt[Cos[c + d*x]]*(a + a*Sec
[c + d*x])^(3/2)) + ((5*A - B)*Sin[c + d*x])/(2*a*d*Sqrt[Cos[c + d*x]]*Sqrt[a + a*Sec[c + d*x]])
Rule 2955
Int[((a_.) + csc[(e_.) + (f_.)*(x_)]*(b_.))^(m_.)*(csc[(e_.) + (f_.)*(x_)]*(d_.) + (c_))^(n_.)*((g_.)*sin[(e_.
) + (f_.)*(x_)])^(p_.), x_Symbol] :> Dist[(g*Csc[e + f*x])^p*(g*Sin[e + f*x])^p, Int[((a + b*Csc[e + f*x])^m*(
c + d*Csc[e + f*x])^n)/(g*Csc[e + f*x])^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, n, p}, x] && NeQ[b*c - a*d
, 0] && !IntegerQ[p] && !(IntegerQ[m] && IntegerQ[n])
Rule 4020
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> -Simp[((A*b - a*B)*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*(d*Csc[e + f*x])^n)/(b*f*(2
*m + 1)), x] - Dist[1/(a^2*(2*m + 1)), Int[(a + b*Csc[e + f*x])^(m + 1)*(d*Csc[e + f*x])^n*Simp[b*B*n - a*A*(2
*m + n + 1) + (A*b - a*B)*(m + n + 1)*Csc[e + f*x], x], x], x] /; FreeQ[{a, b, d, e, f, A, B, n}, x] && NeQ[A*
b - a*B, 0] && EqQ[a^2 - b^2, 0] && LtQ[m, -2^(-1)] && !GtQ[n, 0]
Rule 4013
Int[(csc[(e_.) + (f_.)*(x_)]*(d_.))^(n_)*(csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_))^(m_)*(csc[(e_.) + (f_.)*(x_)]*
(B_.) + (A_)), x_Symbol] :> Simp[(A*Cot[e + f*x]*(a + b*Csc[e + f*x])^m*(d*Csc[e + f*x])^n)/(f*n), x] - Dist[(
a*A*m - b*B*n)/(b*d*n), Int[(a + b*Csc[e + f*x])^m*(d*Csc[e + f*x])^(n + 1), x], x] /; FreeQ[{a, b, d, e, f, A
, B, m, n}, x] && NeQ[A*b - a*B, 0] && EqQ[a^2 - b^2, 0] && EqQ[m + n + 1, 0] && !LeQ[m, -1]
Rule 3808
Int[Sqrt[csc[(e_.) + (f_.)*(x_)]*(d_.)]/Sqrt[csc[(e_.) + (f_.)*(x_)]*(b_.) + (a_)], x_Symbol] :> Dist[(-2*b*d)
/(a*f), Subst[Int[1/(2*b - d*x^2), x], x, (b*Cot[e + f*x])/(Sqrt[a + b*Csc[e + f*x]]*Sqrt[d*Csc[e + f*x]])], x
] /; FreeQ[{a, b, d, e, f}, x] && EqQ[a^2 - b^2, 0]
Rule 206
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1*ArcTanh[(Rt[-b, 2]*x)/Rt[a, 2]])/(Rt[a, 2]*Rt[-b, 2]), x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])
Rubi steps
\begin{align*} \int \frac{\sqrt{\cos (c+d x)} (A+B \sec (c+d x))}{(a+a \sec (c+d x))^{3/2}} \, dx &=\left (\sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{A+B \sec (c+d x)}{\sqrt{\sec (c+d x)} (a+a \sec (c+d x))^{3/2}} \, dx\\ &=-\frac{(A-B) \sin (c+d x)}{2 d \sqrt{\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac{\left (\sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{\frac{1}{2} a (5 A-B)-a (A-B) \sec (c+d x)}{\sqrt{\sec (c+d x)} \sqrt{a+a \sec (c+d x)}} \, dx}{2 a^2}\\ &=-\frac{(A-B) \sin (c+d x)}{2 d \sqrt{\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac{(5 A-B) \sin (c+d x)}{2 a d \sqrt{\cos (c+d x)} \sqrt{a+a \sec (c+d x)}}-\frac{\left ((7 A-3 B) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \int \frac{\sqrt{\sec (c+d x)}}{\sqrt{a+a \sec (c+d x)}} \, dx}{4 a}\\ &=-\frac{(A-B) \sin (c+d x)}{2 d \sqrt{\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac{(5 A-B) \sin (c+d x)}{2 a d \sqrt{\cos (c+d x)} \sqrt{a+a \sec (c+d x)}}+\frac{\left ((7 A-3 B) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}\right ) \operatorname{Subst}\left (\int \frac{1}{2 a-x^2} \, dx,x,-\frac{a \sqrt{\sec (c+d x)} \sin (c+d x)}{\sqrt{a+a \sec (c+d x)}}\right )}{2 a d}\\ &=-\frac{(7 A-3 B) \tanh ^{-1}\left (\frac{\sqrt{a} \sqrt{\sec (c+d x)} \sin (c+d x)}{\sqrt{2} \sqrt{a+a \sec (c+d x)}}\right ) \sqrt{\cos (c+d x)} \sqrt{\sec (c+d x)}}{2 \sqrt{2} a^{3/2} d}-\frac{(A-B) \sin (c+d x)}{2 d \sqrt{\cos (c+d x)} (a+a \sec (c+d x))^{3/2}}+\frac{(5 A-B) \sin (c+d x)}{2 a d \sqrt{\cos (c+d x)} \sqrt{a+a \sec (c+d x)}}\\ \end{align*}
Mathematica [A] time = 1.8762, size = 198, normalized size = 1.12 \[ \frac{2 \tan (c+d x) \sqrt{1-\sec (c+d x)} \left (2 A^2 \cos (2 (c+d x))+2 A^2+A (5 A+3 B) \cos (c+d x)+5 A B-B^2\right )+4 \sqrt{2} (7 A-3 B) \sin \left (\frac{1}{2} (c+d x)\right ) \cos ^3\left (\frac{1}{2} (c+d x)\right ) \sec ^{\frac{3}{2}}(c+d x) \tan ^{-1}\left (\frac{\sqrt{2} \sqrt{\sec (c+d x)}}{\sqrt{1-\sec (c+d x)}}\right ) (A \cos (c+d x)+B)}{4 d \sqrt{\cos (c+d x)-1} (a (\sec (c+d x)+1))^{3/2} (A \cos (c+d x)+B)} \]
Warning: Unable to verify antiderivative.
[In]
Integrate[(Sqrt[Cos[c + d*x]]*(A + B*Sec[c + d*x]))/(a + a*Sec[c + d*x])^(3/2),x]
[Out]
(4*Sqrt[2]*(7*A - 3*B)*ArcTan[(Sqrt[2]*Sqrt[Sec[c + d*x]])/Sqrt[1 - Sec[c + d*x]]]*Cos[(c + d*x)/2]^3*(B + A*C
os[c + d*x])*Sec[c + d*x]^(3/2)*Sin[(c + d*x)/2] + 2*(2*A^2 + 5*A*B - B^2 + A*(5*A + 3*B)*Cos[c + d*x] + 2*A^2
*Cos[2*(c + d*x)])*Sqrt[1 - Sec[c + d*x]]*Tan[c + d*x])/(4*d*Sqrt[-1 + Cos[c + d*x]]*(B + A*Cos[c + d*x])*(a*(
1 + Sec[c + d*x]))^(3/2))
________________________________________________________________________________________
Maple [A] time = 0.385, size = 235, normalized size = 1.3 \begin{align*}{\frac{-1+\cos \left ( dx+c \right ) }{2\,d{a}^{2} \left ( \sin \left ( dx+c \right ) \right ) ^{3}} \left ( 4\,A \left ( \cos \left ( dx+c \right ) \right ) ^{2}\sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}}+A\cos \left ( dx+c \right ) \sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}}+7\,A\sin \left ( dx+c \right ) \arctan \left ( 1/2\,\sin \left ( dx+c \right ) \sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}} \right ) -B\cos \left ( dx+c \right ) \sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}}-3\,B\sin \left ( dx+c \right ) \arctan \left ( 1/2\,\sin \left ( dx+c \right ) \sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}} \right ) -5\,A\sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}}+B\sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}} \right ) \sqrt{\cos \left ( dx+c \right ) }\sqrt{{\frac{a \left ( \cos \left ( dx+c \right ) +1 \right ) }{\cos \left ( dx+c \right ) }}}{\frac{1}{\sqrt{-2\, \left ( \cos \left ( dx+c \right ) +1 \right ) ^{-1}}}}} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
int((A+B*sec(d*x+c))*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x)
[Out]
1/2/d*(-1+cos(d*x+c))*(4*A*cos(d*x+c)^2*(-2/(cos(d*x+c)+1))^(1/2)+A*cos(d*x+c)*(-2/(cos(d*x+c)+1))^(1/2)+7*A*s
in(d*x+c)*arctan(1/2*sin(d*x+c)*(-2/(cos(d*x+c)+1))^(1/2))-B*cos(d*x+c)*(-2/(cos(d*x+c)+1))^(1/2)-3*B*sin(d*x+
c)*arctan(1/2*sin(d*x+c)*(-2/(cos(d*x+c)+1))^(1/2))-5*A*(-2/(cos(d*x+c)+1))^(1/2)+B*(-2/(cos(d*x+c)+1))^(1/2))
*cos(d*x+c)^(1/2)*(a*(cos(d*x+c)+1)/cos(d*x+c))^(1/2)/a^2/(-2/(cos(d*x+c)+1))^(1/2)/sin(d*x+c)^3
________________________________________________________________________________________
Maxima [B] time = 2.42527, size = 11081, normalized size = 62.96 \begin{align*} \text{result too large to display} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((A+B*sec(d*x+c))*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x, algorithm="maxima")
[Out]
-1/4*((4*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*
d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x +
3/2*c)^4 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^4 + 4*(7*log(cos(
1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin
(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^4 + 70*(log(c
os(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + si
n(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c)^2 + 7*(log(cos
(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(
1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^4 - 8*sin(1/2*d*x + 1/2*c)^5 + 28*(7*(l
og(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2
+ sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*
d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^3 + 4*(21*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*
d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2
*d*x + 1/2*c) - 24*sin(1/2*d*x + 1/2*c)^2 - 20)*sin(3/2*d*x + 3/2*c)^3 - 8*(10*cos(1/2*d*x + 1/2*c)^2 + 3)*sin
(1/2*d*x + 1/2*c)^3 + ((7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) -
7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*
cos(3/2*d*x + 3/2*c)^2 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1)
- log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2 +
(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/
2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2
+ 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/
2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 8*sin(1/2*d*x + 1/2*c)
^3 + 6*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x
+ 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2
*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) + 2*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2
*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1)
)*sin(1/2*d*x + 1/2*c) - 8*sin(1/2*d*x + 1/2*c)^2 - 8)*sin(3/2*d*x + 3/2*c) - 8*(9*cos(1/2*d*x + 1/2*c)^2 + 2)
*sin(1/2*d*x + 1/2*c))*cos(5/2*d*x + 5/2*c)^2 + (427*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*
sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))
*cos(1/2*d*x + 1/2*c)^2 + 35*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1
) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2
- 40*sin(1/2*d*x + 1/2*c)^3 - 8*(61*cos(1/2*d*x + 1/2*c)^2 + 9)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^2 +
((7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1
/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^
2 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x +
1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2 + (7*log(cos(1/2*d*x +
1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x
+ 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2 + 7*(log(cos(1/2*d*x
+ 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x
+ 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 8*sin(1/2*d*x + 1/2*c)^3 + 6*(7*(log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2
*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2
*c))*cos(3/2*d*x + 3/2*c) + 2*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c)
+ 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c
) - 8*sin(1/2*d*x + 1/2*c)^2 - 8)*sin(3/2*d*x + 3/2*c) - 8*(9*cos(1/2*d*x + 1/2*c)^2 + 2)*sin(1/2*d*x + 1/2*c)
)*sin(5/2*d*x + 5/2*c)^2 + (8*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c)
+ 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1
/2*c))*cos(3/2*d*x + 3/2*c)^2 + 259*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2
*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/
2*c)^2 + 91*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d
*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 104*sin(1/2*d*x
+ 1/2*c)^3 + 28*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(c
os(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2
*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) - 8*(37*cos(1/2*d*x + 1/2*c)^2 + 21)*sin(1/2*d*x + 1/
2*c))*sin(3/2*d*x + 3/2*c)^2 + 2*(2*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1
/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d
*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^3 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x
+ 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x
+ 1/2*c)^3 + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1
/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1
/2*c)^2 - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^3 + 13*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/
2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1
/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)^2 + (2*(7
*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c
)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) + 7*
(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^
2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/
2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2 + 2*(84*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/
2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1
/2*d*x + 1/2*c)^2 + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log
(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 8*sin
(1/2*d*x + 1/2*c)^3 - 16*(6*cos(1/2*d*x + 1/2*c)^2 + 1)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) + 2*(7*(log
(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 +
sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c) - 8*cos(1/2*d*
x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 + 2*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x
+ 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x
+ 1/2*c) - 8*sin(1/2*d*x + 1/2*c)^2 - 8)*cos(3/2*d*x + 3/2*c) - 8*cos(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c)
- 8*(9*cos(1/2*d*x + 1/2*c)^3 + 2*cos(1/2*d*x + 1/2*c))*sin(1/2*d*x + 1/2*c))*cos(5/2*d*x + 5/2*c) + 2*(147*(l
og(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2
+ sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^3 + 35*(log(cos(1/2*d*x + 1/2*c)^
2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2
- 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 - 40*cos(1/2*d*x + 1/2*c)*sin(1/2*
d*x + 1/2*c)^3 - 56*(3*cos(1/2*d*x + 1/2*c)^3 + cos(1/2*d*x + 1/2*c))*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*
c) + 2*(2*(7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 7*log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x +
3/2*c)^3 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x +
1/2*c) + 7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*
x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^3 - 8*sin(1/2*d*x +
1/2*c)^4 + (7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2
*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c) - 8*sin(1/2*d*x +
1/2*c)^2 - 4)*cos(3/2*d*x + 3/2*c)^2 + (35*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d
*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*
d*x + 1/2*c) - 40*sin(1/2*d*x + 1/2*c)^2 - 36)*sin(3/2*d*x + 3/2*c)^2 - 4*(18*cos(1/2*d*x + 1/2*c)^2 + 5)*sin(
1/2*d*x + 1/2*c)^2 + 6*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) -
log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c)*sin(1
/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 - 4*cos(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c)
- 36*cos(1/2*d*x + 1/2*c)^2 + 2*((7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2
*c) + 1) - 7*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 8*sin(1/2*d*x
+ 1/2*c))*cos(3/2*d*x + 3/2*c)^2 + 63*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x +
1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x +
1/2*c)^2 + 14*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/
2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(1/2*d*x + 1/2*c)^2 - 16*sin(1/2*d
*x + 1/2*c)^3 + 6*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(
cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2*d*x + 1/2*c) - 8*cos(1/
2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*cos(3/2*d*x + 3/2*c) - 4*(18*cos(1/2*d*x + 1/2*c)^2 + 7)*sin(1/2*d*x + 1/
2*c))*sin(3/2*d*x + 3/2*c))*sin(5/2*d*x + 5/2*c) + 2*(133*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2
+ 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c)
+ 1))*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c) + 21*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 +
2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1
))*sin(1/2*d*x + 1/2*c)^3 - 24*sin(1/2*d*x + 1/2*c)^4 + 2*(21*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*
c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2
*c) + 1))*sin(1/2*d*x + 1/2*c) - 24*sin(1/2*d*x + 1/2*c)^2 - 20)*cos(3/2*d*x + 3/2*c)^2 - 8*(19*cos(1/2*d*x +
1/2*c)^2 + 7)*sin(1/2*d*x + 1/2*c)^2 + 16*(7*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*
d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(1/2
*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c) - 8*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 - 5*cos(1/2*d*x + 1/2*c))*c
os(3/2*d*x + 3/2*c) - 80*cos(1/2*d*x + 1/2*c)^2)*sin(3/2*d*x + 3/2*c) - 8*(9*cos(1/2*d*x + 1/2*c)^4 + 11*cos(1
/2*d*x + 1/2*c)^2)*sin(1/2*d*x + 1/2*c))*A*sqrt(a)/(4*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^4 + 28*sqrt(2)*a^2*cos(
3/2*d*x + 3/2*c)^3*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^4 + 4*sqrt(2)*a^2*sin(3/2*d*x + 3
/2*c)^4 + 12*sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)^3*sin(1/2*d*x + 1/2*c) + 10*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2*s
in(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^4 + (sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2 + 6*sqrt(2)*a
^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(3/2*d*x
+ 3/2*c)^2 + 2*sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*cos
(5/2*d*x + 5/2*c)^2 + (61*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + 5*sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*cos(3/2*d
*x + 3/2*c)^2 + (sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2 + 6*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c)
+ 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)^2 + 2*sqrt(2)*a^2*sin(3/2*d*x + 3/2*
c)*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*sin(5/2*d*x + 5/2*c)^2 + (8*sqrt(2)*a^2*cos(3/2*
d*x + 3/2*c)^2 + 28*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c) + 37*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c
)^2 + 13*sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*sin(3/2*d*x + 3/2*c)^2 + 2*(2*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^3
+ 13*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^3 + sqrt(2)*
a^2*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2 + (2*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c) + sqrt(2)*a^2*cos(1/2*d*
x + 1/2*c))*sin(3/2*d*x + 3/2*c)^2 + 2*(12*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2 + sqrt(2)*a^2*sin(1/2*d*x + 1/2*
c)^2)*cos(3/2*d*x + 3/2*c) + 2*(2*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*cos(1/2*
d*x + 1/2*c)*sin(1/2*d*x + 1/2*c))*sin(3/2*d*x + 3/2*c))*cos(5/2*d*x + 5/2*c) + 2*(21*sqrt(2)*a^2*cos(1/2*d*x
+ 1/2*c)^3 + 5*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c)^2)*cos(3/2*d*x + 3/2*c) + 2*(2*sqrt(2)*a^
2*sin(3/2*d*x + 3/2*c)^3 + sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)^2*sin(1/2*d*x + 1/2*c) + 6*sqrt(2)*a^2*cos(3/2*d*x
+ 3/2*c)*cos(1/2*d*x + 1/2*c)*sin(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c
) + 5*sqrt(2)*a^2*sin(3/2*d*x + 3/2*c)^2*sin(1/2*d*x + 1/2*c) + sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^3 + 2*(sqrt(2
)*a^2*cos(3/2*d*x + 3/2*c)^2 + 6*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c) + 9*sqrt(2)*a^2*cos(1/2
*d*x + 1/2*c)^2 + 2*sqrt(2)*a^2*sin(1/2*d*x + 1/2*c)^2)*sin(3/2*d*x + 3/2*c))*sin(5/2*d*x + 5/2*c) + 2*(6*sqrt
(2)*a^2*cos(3/2*d*x + 3/2*c)^2*sin(1/2*d*x + 1/2*c) + 16*sqrt(2)*a^2*cos(3/2*d*x + 3/2*c)*cos(1/2*d*x + 1/2*c)
*sin(1/2*d*x + 1/2*c) + 19*sqrt(2)*a^2*cos(1/2*d*x + 1/2*c)^2*sin(1/2*d*x + 1/2*c) + 3*sqrt(2)*a^2*sin(1/2*d*x
+ 1/2*c)^3)*sin(3/2*d*x + 3/2*c)) - (3*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x +
1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(2*d*x +
2*c)^2 + 12*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d
*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(d*x + c)^2 + 3*(log(cos(1/2*d*x + 1/
2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2
*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(2*d*x + 2*c)^2 + 12*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c
)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*
c) + 1))*sin(d*x + c)^2 + 2*(6*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) +
1) - log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*cos(d*x + c) + 3*log(
cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 3*log(cos(1/2*d*x + 1/2*c)^2 +
sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) - 2*sin(3/2*d*x + 3/2*c) + 2*sin(1/2*d*x + 1/2*c))*cos(2
*d*x + 2*c) + 4*(3*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - 3*log(c
os(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1) + 2*sin(1/2*d*x + 1/2*c))*cos(d*x
+ c) + 4*(3*(log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*sin(1/2*d*x + 1/2*c) + 1) - log(cos(1/2*
d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1))*sin(d*x + c) + cos(3/2*d*x + 3/2*c) - c
os(1/2*d*x + 1/2*c))*sin(2*d*x + 2*c) - 4*(2*cos(d*x + c) + 1)*sin(3/2*d*x + 3/2*c) + 8*cos(3/2*d*x + 3/2*c)*s
in(d*x + c) - 8*cos(1/2*d*x + 1/2*c)*sin(d*x + c) + 3*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 + 2*
sin(1/2*d*x + 1/2*c) + 1) - 3*log(cos(1/2*d*x + 1/2*c)^2 + sin(1/2*d*x + 1/2*c)^2 - 2*sin(1/2*d*x + 1/2*c) + 1
) + 4*sin(1/2*d*x + 1/2*c))*B/((sqrt(2)*a*cos(2*d*x + 2*c)^2 + 4*sqrt(2)*a*cos(d*x + c)^2 + sqrt(2)*a*sin(2*d*
x + 2*c)^2 + 4*sqrt(2)*a*sin(2*d*x + 2*c)*sin(d*x + c) + 4*sqrt(2)*a*sin(d*x + c)^2 + 4*sqrt(2)*a*cos(d*x + c)
+ 2*(2*sqrt(2)*a*cos(d*x + c) + sqrt(2)*a)*cos(2*d*x + 2*c) + sqrt(2)*a)*sqrt(a)))/d
________________________________________________________________________________________
Fricas [A] time = 0.53, size = 1068, normalized size = 6.07 \begin{align*} \left [-\frac{\sqrt{2}{\left ({\left (7 \, A - 3 \, B\right )} \cos \left (d x + c\right )^{2} + 2 \,{\left (7 \, A - 3 \, B\right )} \cos \left (d x + c\right ) + 7 \, A - 3 \, B\right )} \sqrt{a} \log \left (-\frac{a \cos \left (d x + c\right )^{2} - 2 \, \sqrt{2} \sqrt{a} \sqrt{\frac{a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt{\cos \left (d x + c\right )} \sin \left (d x + c\right ) - 2 \, a \cos \left (d x + c\right ) - 3 \, a}{\cos \left (d x + c\right )^{2} + 2 \, \cos \left (d x + c\right ) + 1}\right ) - 4 \,{\left (4 \, A \cos \left (d x + c\right ) + 5 \, A - B\right )} \sqrt{\frac{a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt{\cos \left (d x + c\right )} \sin \left (d x + c\right )}{8 \,{\left (a^{2} d \cos \left (d x + c\right )^{2} + 2 \, a^{2} d \cos \left (d x + c\right ) + a^{2} d\right )}}, \frac{\sqrt{2}{\left ({\left (7 \, A - 3 \, B\right )} \cos \left (d x + c\right )^{2} + 2 \,{\left (7 \, A - 3 \, B\right )} \cos \left (d x + c\right ) + 7 \, A - 3 \, B\right )} \sqrt{-a} \arctan \left (\frac{\sqrt{2} \sqrt{-a} \sqrt{\frac{a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt{\cos \left (d x + c\right )}}{a \sin \left (d x + c\right )}\right ) + 2 \,{\left (4 \, A \cos \left (d x + c\right ) + 5 \, A - B\right )} \sqrt{\frac{a \cos \left (d x + c\right ) + a}{\cos \left (d x + c\right )}} \sqrt{\cos \left (d x + c\right )} \sin \left (d x + c\right )}{4 \,{\left (a^{2} d \cos \left (d x + c\right )^{2} + 2 \, a^{2} d \cos \left (d x + c\right ) + a^{2} d\right )}}\right ] \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((A+B*sec(d*x+c))*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x, algorithm="fricas")
[Out]
[-1/8*(sqrt(2)*((7*A - 3*B)*cos(d*x + c)^2 + 2*(7*A - 3*B)*cos(d*x + c) + 7*A - 3*B)*sqrt(a)*log(-(a*cos(d*x +
c)^2 - 2*sqrt(2)*sqrt(a)*sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))*sin(d*x + c) - 2*a*cos(d*
x + c) - 3*a)/(cos(d*x + c)^2 + 2*cos(d*x + c) + 1)) - 4*(4*A*cos(d*x + c) + 5*A - B)*sqrt((a*cos(d*x + c) + a
)/cos(d*x + c))*sqrt(cos(d*x + c))*sin(d*x + c))/(a^2*d*cos(d*x + c)^2 + 2*a^2*d*cos(d*x + c) + a^2*d), 1/4*(s
qrt(2)*((7*A - 3*B)*cos(d*x + c)^2 + 2*(7*A - 3*B)*cos(d*x + c) + 7*A - 3*B)*sqrt(-a)*arctan(sqrt(2)*sqrt(-a)*
sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))/(a*sin(d*x + c))) + 2*(4*A*cos(d*x + c) + 5*A - B)*
sqrt((a*cos(d*x + c) + a)/cos(d*x + c))*sqrt(cos(d*x + c))*sin(d*x + c))/(a^2*d*cos(d*x + c)^2 + 2*a^2*d*cos(d
*x + c) + a^2*d)]
________________________________________________________________________________________
Sympy [F(-1)] time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((A+B*sec(d*x+c))*cos(d*x+c)**(1/2)/(a+a*sec(d*x+c))**(3/2),x)
[Out]
Timed out
________________________________________________________________________________________
Giac [F] time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{{\left (B \sec \left (d x + c\right ) + A\right )} \sqrt{\cos \left (d x + c\right )}}{{\left (a \sec \left (d x + c\right ) + a\right )}^{\frac{3}{2}}}\,{d x} \end{align*}
Verification of antiderivative is not currently implemented for this CAS.
[In]
integrate((A+B*sec(d*x+c))*cos(d*x+c)^(1/2)/(a+a*sec(d*x+c))^(3/2),x, algorithm="giac")
[Out]
integrate((B*sec(d*x + c) + A)*sqrt(cos(d*x + c))/(a*sec(d*x + c) + a)^(3/2), x)